NASA 

Technical 

Paper 

1999 

May 1982 


NASA 

Technical 

Paper 

1999 


1982 


TECH LIBRARY KAFB, NM 


iiiniiiilii 

□□taibi 


Computer Modeling 
of Multiple-Channel Input 
Signals and Intermodulation 
Losses Caused by Nonlinear 
Traveling-Wave-Tube Amplifiers 


N. Stankiewicz 
Lewis Research Center 
Cleveland, Ohio 


rviASA 

National Aeronautics 
and Space Administration 


Scientific and Technical 
Information Branch 


Summary 

The multiple-channel input signal to a soft limiter 
amplifier such as a traveling wave tube is represented as a 
finite, linear sum of Gaussian functions in the frequency 
domain. Linear regression is used to fit the channel 
shapes to a least-squares residual error. 

Distortions in output signal, namely intermodulation 
products, are produced by the nonlinear gain 
characteristic of the amplifier and constitute the principle 
noise analyzed in this study. The signal-to-noise ratios are 
calculated for various input powers from saturation to 10 
dB below saturation for two specific distributions of 
channels. 

A criterion for the truncation of the series expansion of 
the nonlinear transfer characteristic is given. It is found 
that the signal-to-noise ratios are very sensitive to the 
coefficients used in this expansion. Improper or incorrect 
truncation of the series leads to ambiguous results in the 
signal-to-noise ratios. 

Introduction 

The multichannel use of an amplifier in a 
communications system requires linear operation in order 
to prevent interference between channels. For a soft- 
limiting, nonlinear device such as a traveling wave tube 
(TWT) amplifier this means backing off below saturation 
to the linear portion of the gain curve. Such quasi-linear 
operation results in an unavoidable loss of efficiency 
which can be partially recovered in a TWT through the 
use of a depressed collector. 

Despite the advances (ref. 1) of multistaged depressed 
collectors in recovering the spent beam energy of TWT’s, 
it is important to determine the extent of the backoff if 
the highest possible efficiency is to be achieved. This 
importance manifests itself in two ways. First, the 
collector design is a function of the TWT operating 
point. And second, the TWT gain characteristics may be 
a design option to be traded with other parameters; for 
exEunple, the length of the slow wave structure. 

The principal nonlinearity of TWT’s affecting 
efficiency is the power gain characteristic (AM/AM 
conversion). This type of nonlinearity wastes power 
directly by producing intermodulation bands and by 
broadening the channel signals. These out-of-band 
signals must be removed with filters and represent a loss 
in available power. 

Other nonlinearities such as amplitude-to-phase 
conversion (AM/PM) also affect efficiency but indirectly 
through distortion of the input signal. Acceptable signal 
quality depends upon the type of modulation and 
demodulation to be used. Such considerations are beyond 
the scope of this paper. It will be assumed that operating 
quasi-linearly to diminish AM/AM conversion will also 


diminish AM/PM distortions to an acceptable level. 
Therefore, only the AM/AM characteristics will be 
discussed in this report. 

The purpose of the work reported in this paper is to 
develop a computational program to evaluate AM/AM 
channel interference in a TWT. The computer program is 
designed to use empirical data in the modeling of the 
input channel spectra and also for the gain 
characteristics. 

By channel spectra we mean the amplitude distribution 
of the channels as a function of frequency. Early in the 
development of this computer program it was decided 
that modeling of input spectrum was a desirable option 
for several reasons. First, this approach provides 
flexibility in analyzing input signals without direct 
consideration of the method of modulation or signal 
coding that may be used in the communications system. 
Second, questions concerning accurate representations of 
signals which convey real information and which are 
usually handled by taking ensemble averages (ref. 2) over 
these randomly occurring variables can be circumvented. 
Finally, the method chosen for the modeling of the 
channels, a least-squares expansion in Gaussian 
functions, is reminiscent of the central limit theorem in 
which the cumulative effect of many random processes 
tends to a Gaussian distribution. 

It was envisioned that this work would be 
supplemented by an auxiliary systems program which 
would simulate a modulation technique, perform an 
ensemble average over a random signal (both 
information as well as noise) and spectrally analyze the 
result. This or, perhaps, measured spectral data would be 
the input for the AM/ AM interference program. 

The voltage gain curve must also be modeled 
realistically from measured or computed data. 
Mathematically, it is an odd function and as such can be 
expanded in odd powers of the input voltage. The 
coefficients of the polynomial fit are calculated from a 
least-squares algorithm. The deviation from linearity 
arises, of course, from the thirds, fifths, etc., terms of the 
expression; and each additional term gives rise to a higher 
order intermodulation product. A problem arises in 
deciding where to truncate the series. From the point of 
view of the least-squares error the fit becomes better as 
higher order terms are included. However, with each 
additional term the expansion coefficients must be 
recalculated and will differ from previously calculated 
coefficients thus causing an ambiguity in the 
computation of intermodulation products. This is a 
consequence of the nonorthogonality of the power series. 
A method for removing the ambiguity consistent with the 
inherent error of the data is discussed in appendix A. 

This paper is divided into a section called 
FORMULATION and one called RESULTS with most 
of the mathematical derivations relegated to appendixes 



A, B, and C. The computer program is not given at this 
time but is demonstrated in the section RESULTS using 
the particular channel spacings of reference 3. 

\ 

Formulation 

Figure 1 is typical of a power gain curve for a TWT. 
The input and output powers have been normalized to 
their saturation values and the deviation from linearity 
(dashed line) is the gain compression (about 1.9 dB at 
saturation). This particular curve is the result of TWT 
computer simulation and was generated as part of a 
design study for the 75-watt, 20-gigahertz tube of 
reference 3. 

The voltage transfer curve which is the square root of 
the power gain curve is shown in figure 2. The error 
associated with the data used in making this plot was 
judged to be ±0.005 based on visually estimating the 
points in figure 1 . The error would have been smaller had 
the raw computer data been used. However, experimental 
data is usually plotted and the experimental accuracy 
often is not given; hence, the error associated with a 
visual estimate is probably consistent with most 
generated data. 

The method described in appendix A can be used to 
expand the vector of voltage transfer data in an 
orthogonal set of polynomials fVj 

(^/^«t)output = r^y»0- 0) 



Figure 1. - Normalized power gain curve for a 
traveling wave tube. 



Figure 2. - Normalized voltage transfer curve for a travel- 
ing wave tube. 


A computer program called POLYBASE was written to 
compute the orthogonal vectors fVj from the input vector 
data (^^/f^sat)input perform the scalar products 

that define the expansion coefficients bj as follows: 

= ( 2 ) 

The results obtained by using the data of figure 2 are 
given in table I for the first six coefficients. Only the first 
three coefficients fall within the error norm of 
2.022 X 10 “2 as calculated from this data using equation 
(A8). Therefore, the data should be fit with a fifth-order 
polynomial. By using a least-squares-error algorithm, the 
fifth-order polynomial fit of the data in figure 2 is found 
to be 


j=1.228x + 0.158x-3 -0.386x5 (3) 

TABLE 1. - FIRST SIX SCALAR PRODUCTS OF 
OUTPUT DATA VECTOR (SEE FIG. 2) WITH 
ORTHOGONAL BASIS OF DATA 
VECTOR SPACE 


*1 

3.145 

bi 

-0.2430 

b. 

-0.8118x10-' 


— Error norm = 0.2022 x 10 “ ' 

b. 

-0.3194x10“^ 

bs 

0.9272 X 10"^ 

bt 

0.1286x10“^ 
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where y = ( K/ and x= ( V / . 

It will be assumed that equation (3) is valid when the 
input and output signals are functions of time. In 
general, the output signal appears at a time t after 
introducing the input signal. That is, an output y(t) will 
occur in response to an input xU + t). However, unless 
the time delay t is a function of frequency, the Fourier 
transform of equation (3) will yield a constant and hence 
factorable phase for each term in the equation. The study 
of frequency dependent time delays is an interesting facet 
of nonlinear systems; however, as explained in the 
INTRODUCTION, they are beyond the scope of this 
paper. 

Equation (3) can thus be rewritten as 

y(t) = 2] C 2 *_ ' (4) 

k = ! 



Figure 3. - Least-squares fit of a trapezoidai frequency gate 
using seven Gaussian functions. 


The Fourier transform of equation (4) gives the following 
equivalent expression in the frequency domain (see eq. 
(B4)): 




(5) 


The normalized output spectral power density from 
equation (B42) is 


>^C/)= t t C2k^-XC2k,-X 

A:i = 1 k2~ \ 


[^(/)®(2A. - D^(y)®(2A:, - 1) 


where 

P(k,.k^) = - D^,3tC/)®(2A:2- Drf/ (6) 


The input signal X(f) to the nonlinear system is assumed 
to be a distributed rather than a line source. The 
distribution can be expanded in a series of Gaussitm 
sources each having a width, that is, a standard deviation 
of 1/V2X where X is a convenient parameter to adjust the 
expansion. Equation (B13) or its integerized version (eq. 
(C24)) is the type of expansion that is supposed. Using 
only the positive frequencies from equation (B13) we 
have 


X(f)^ 





(7) 


Equation (7) is a linear expansion in N Gaussian 
functions each centered about a frequency f„ with the 
amplitudes a„ as the expansion coefficients. A linear 
regression with a least -squares technique can be used to 
find a„ for a given distribution. Figure 3 shows the results 
of such a computation for a trapezoidal frequency gate or 
bandpass filter. The gate is ideal and has unit amplitude 
at the discrete relative frequencies between 5 to 9 and zero 
amplitude elsewhere. The calculation was done for two 
values of the parameter X. The solid curve with X= 1 has 
the smaller residual and therefore is chosen as the model 
of a bandpass filter. In the remainder of this paper the 
input driving signal to the nonlinear amplifier will be in 
the form of multiple channels defined by bandpass filter 
gates. 

The amplitude coefficients have, so far, been 
arbitrary and it is necessary to normalize them in order to 
compute the correct results. The input power P; is 
defined by the integral 


p,=r x{f)^df 

J — oo 


( 8 ) 


Using equation (7) we get 
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After rearranging terms and carrying out the integration, 
this becomes 





-(X/2)(/„, 


rt|=l rt2=l 


( 10 ) 


The scale factor 5 from equation (B39) can be used to 
rewrite equation (10) as 




-(X/2)0;, -/„,)= 


Results 

It was assumed in reference 3 that the total bandwidth 
available to a certain communications system was divided 
into nine separate channels of equal width (see fig. 4). 
The channels are delineated with ideal bandpass filters. It 
was further assumed that the transmitter would consist of 
three TWT amplifiers carrying three channnels each. The 
distribution of the three channels within each TWT was 
chosen so that the third-order intermodulation products 
would fall between the channels and could, therefore, be 
removed with filters after amplification. 

Assigning channel numbers (as in fig. 4) means, the 
first TWT would carry channels 1,3, and 7; the second 
TWT would carry channels 2, 4, and 8; and the third 



( 11 ) 


Equation (1 1) becomes the defining equation for the scale 
factor S, which is equal to unity at saturation and less 
than one below saturation. The arbitrariness in the 
amplitude coefficients a„ is removed by redefining them 
as 


an-a„/{PiYir (12) 

A computer program was developed to calculate the 
output spectral density (eq. 5)) as a function of the input 
signal (eq. (7)). The program coding, which will not be 
presented, follows the logic outlined in appendix C. The 
coherent part of the spectral density T^ is computed from 
the terms given by equations (C24) to (C29). The noise 
density is then given by 

Y%= _ y2 (13) 
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Figure 4. - Assumed distribution of nine channeis over 
three traveiing wave tubes (TWT's). 
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Figure 5. - Distribution of channels for case one. 


The signal-to-noise ratio in the band between/^ and/* is 

(5/A)„,* = j^*y2(/)// J^V^(/)c/ (14) 

As an example, the program was used to aneilyze the 
signal-to-noise ratios of the multichannel inputs to a 
nonlinear TWT as proposed in reference 3. These 
computations are discussed in the following section. 



Figure 6. - Distribution of channeis for case two. 
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Normalized output power Normalized input power Normalized output power Normalized input power 




TWT would carry the remaining channels 5, 6, and 9. If it 
is assumed that the bandwidth of each TWT is uniform, 
then the channel distribution 1,3, and 7 is identical to the 
distribution 2, 4, and 8. This is because the two 
arrangements are merely translations of each other and 
intermodulation products depend only on relative 


frequencies. Hence, in order to analyze this system only 
two distinct cases need be considered. 

Figure 5 shows the distribution of channels, designated 
as case one, which is the representation in Gaussian 
functions of both of the arrangements (1, 3, 7) and (2, 4, 
8). Figure 6 is the distribution for case two which is the 



J =L 




- 4 ^ 


60 80 100 120 -80 -60 -40 -20 0 20 40 60 80 


TooTlio 


(a) Step: 0 dB (saturation). (b) Step-. -0.5dB. 



(c)Step; -1.0 dB. (d)Step: -l.SdB. 
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(e)Step: -2 0dB. 


(f)Step: -iSdB. 


Figure 11. - Sequence of output noise power for channel distribution of case one in 1/2-dB steps from saturation to 
-5dB below saturation. Power is normaiized to its peak for each step. 
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inputTepresentation for channels 5, 6, and 9. The three 
channels in each case are designated as channels A, B, 
and C. 

Cases one and two were analyzed for a sequence of 
input powers ranging from saturation to - 10 dB below 
saturation. The computer program uses a graphics 


package to display the input and the output spectral 
power density of the nonlinear amplifier. Figure 7 shows 
the power density of the input and the output at 
saturation for case one. As the drive power is decreased 
beyond -3 dB from saturation, the output 
intermodulation structure, although it is present, is 
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(i) Step: -4.0dB. (j) Step: -4.5dB. 



Figure 11. - Concluded. 
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scarcely discernible in the graphics display. The situation 
at - 4 dB is given in figure 8 and shows that the output 
spectral density appears to be a replica of the input 
spectral density. 

Figures 9 and 10 show the spectral density input and 
output power for case two. Figure 9 is the result at 
saturation and figure 10 is at ~4 dB from saturation. 


It is more interesting to have the computer produce a 
display of the noise spectral density rather than the total 
output spectral density. Figure II is a sequence of 
computer plots showing the noise density for case one in 
steps of 0.5 dB from saturation to -5 dB below 
saturation. Figure 12 is an identical sequence but for case 
two. Note that, between -3.5 and -4.5 dB below 



(a) step: OdB (saturation). (b)Step: -0.5dB. 



(c)Step: -LOdB. (d)Step: -1.5dB. 



01 I I /'I I I I I I I I 7i' I '1^1 J'l 'Vi I I 

-30 -20 -10 0 10 20 30 40 50 60 ?0 -30 -20 -10 0 10 20 30 40 50 60 


Relative frequency 

(e)Step: -iOdB. (f) Step: -Z5dB. 

Figure 12. - Sequence of output noise power for channel distribution of case two in 1/2-dB steps from saturation to 
-5 dB below saturation. Power is normalized to its peak for each step. 
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saturation, a type of transition takes place in which the 
character of the noise changes. This behavior is a 
consequence of the particular voltage transfer function 
that was used in this example (see fig. 2). Because the 
least-squares fit to this curve (eq. (3)) has a positive third- 
order coefficient and a negative fifth-order coefficient 


there exists a region in which these two terms cancel. This 
occurs at 


= = 0-409= -3.883 dB 

U.3oO 


(15) 




(k)Slep: -S.OdB. 
Figure 12. - Concluded. 
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Note that at saturation when 1 the fifth-order term 
is larger than the third-order term. Thus, the fifth-order 
terms dominate from saturation to the transition point 
(eq. (15)) and the third-order terms are dominant below 
this point. 

The computer program provides the signal-to-noise 
ratio for each channel by calculating the integrals in 
equation (14). Figures 13 and 14 show the results of this 
computation for the channel distributions of cases one 
and two, respectively. 

The most interesting feature of these plots is the 
apparent singularity at approximately the value given by 
equation (15). This was unexpected because equation (15) 
is the result of integrating over all frequencies while the 
signal-to-noise ratios involve local integrals over the band 
in which a channel is defined. 

The munerical differences in signal-to-noise ratios 
between individual channels and even between the two 
cases are very slight. In figures 13 and 14 the line 
thicknesses bound the values for the various channels. If 
the two cases are superimposed, they coincide within 
these bounds. Evidently the signal-to-noise ratio is very 
dependent on the coefficients of the third- and fifth-order 
terms and the fact that in this case they differ in sign. 

A hypothetical transfer curve in which both of these 
coefficients are negative was derived to study the same 



Normalized input power, dB 


Figure 13. - Signal-to-noise ratio as a func- 
tion of drive power for case one. (Individuai 
channels are bounded by line thicknesses. ) 



Figure 14 - Signal-to-noise ratio as a func- 
tion of drive power for case two. (Individuai 
channels are bounded by line thicknesses. ) 



Figure 15. - Plot of third and fifth coefficients 
as functions of parameter €. 
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channel distributions. The dashed lines in figures 13 and 
14 show these results. 

The voltage transfer curve shown in figure 2 has the 
general characteristics that 


y=i 


Cl =1.3333 
C 3 =0.1667 > 

C 5 =0.1667 

J 


(20) 


y atx=±l 




and 

dy 


J 


dx 


> 1 at j :=0 


(16) 


The first condition is a consequence of using a 
normalized transfer function while the second is the 
saturation characteristic. The third condition is a 
statement that for curves of single curvature one can 
never satisfy conditions ( 1 ) and ( 2 ) unless the coefficient 
of the linear term is greater than one. Using equation (4) 
and the conditions in equation (16) yields 


Cl +C 3 +C 5 = 1 

Ci+C3+C5=0 > 


(17) 


Cl = 1 +€, £>0 J 

These can be solved in terms of e with the following 
result: 


^=-(5) J 


(18) 


Obviously C 3 is always negative for e greater than 1/4 and 
C 5 is always positive for e greater than 1/2. However, 
both coefficients may be negative in the region 

J<‘<5 OS) 


However, the results shown in figures 13 and 14 indicate 
that the hypothetical case is noisier. Further 
investigations of the effect of the coefficients on the 
signal-to-noise ratio were deferred until a more complete 
study could be undertaken. Such an investigation would 
be of great help in designing TWT’s with optimum low 
noise characteristics. 


Concluding Remarks 

The nonlinear power gain characteristic of a soft 
limiter such as a traveling wave tube (TWT) produces 
interchannel distortions if it is used as a multichaimel 
amplifier. In order to minimize the effect of these 
distortions, it is necessary to operate the TWT quasi- 
linearly; that is, backed-off from saturation. With the 
help of a computer program which was developed for this 
problem, the signal-to-noise ratios are calculated for 
various distributions of channels as a function of driving 
power. 

The nonlinear gain characteristic is represented as a 
power series in input signal voltage and a least-squares 
algorithm is used to fit the measured or computed gain 
data. A criterion for truncation of the series is developed 
which depends on the inherent error of the data. This 
criterion removes the ambiguities which results if an 
arbitrary truncation is used. 

The input signal channels are modeled by linear 
regression theory in which measured or computed 
bandpass channel characteristics are fit with Gaussian 
functions. This procedure bypasses the ensemble 
averaging that is usually performed to specify the 
frequency bands required by a given channel. This also 
defers the need to consider signal coding, modulation and 
information recovery in the analysis. 

Two specific distributions of channels are analyzed. It 
is found that the signal-to-noise ratios are very sensitive 
to the coefficients which are used in the power series 
expansion of the gain characteristic curve. The 
procedures described in this paper can be used in defining 
the optimum gain characteristics and drive levels for low 
noise, multichannel operation. 


This region is shown in figure 15 where equations (18) 
are plotted. The values chosen for the hypothetical Lewis Research Center 

transfer function were National Aeronautics and Space Administration 

Cleveland, Ohio, May 28, 1981. 


11 




Appendix A 

Voltage Transfer Function (AM/AM) for a TWT 


We assume that a measured or calculated power 
transfer curve is available. Figure 1 is typical of the 
nonlinear characteristics of such a calculated power 
transfer curve. The normalized output power is plotted in 
this figure as a function of the input or drive power. The 
normalization is with respect to saturation values. 

The normalized voltage transfer function is derived by 
taking the square root of the power curve. This function 
is plotted in figure 2. Note that the voltage transfer 
function includes both the positive and negative roots of 
the power curve and that it is a completely odd function. 

Because it is an odd function, the voltage transfer can 
be represented as a polynomial expansion using only odd 
terms. Defining x(0 zxidyif) as the normalized input and 
output signal voltages yields 

X0= (Al) 

*=i 

The limit on the summation K is the order of the 
truncation of the representation. For example, K=5 will 
truncate the series at the ninth power of the input and this 
will, of course, result in ninth-order intermodulation 
products. 

One way of determining the coefficients (cj, C 3 , C 5 , 

. . . , C 2 k-\) is to apply, directly, a least-squares 
polynomial fit to equation (Al). Let the data, taken from 
the voltage transfer curve, be designated as 


for the coefficients c/. Unfortunately, this method of 
determining the coefficients tells nothing about the 
number of terms required for the representation. That is, 
the final degree of the polynomial is arbitrary. 

The computed deviations from the data show that the 
fit becomes better as more terms are considered. 
However, the coefficients change with the addition of 
each new term. This leads to ambiguous results in which 
each representation gives a different value for the 
intermodulation products. In some radical cases the 
higher order products are more pronounced than the 
lower orders. This cannot be an adequate description of a 
real nonlinear amplifier. 

In order to set a limit to the number of required terms, 
some account must be taken of the quality of data used. 
To begin, equation (Al) can be written as the following 
vector expansion in polynomial space: 


A A 

I 



/vA 

/ n 

1 ^2 


1 xj 

• 

= ci 

\ 

+C 3 

• 

- 1 -C 5 

4- 

w 


w 


W 


W 


yi=y{Xi) (A2) 

where i ranges over the number of data points used in the 
fit. 

The sum of the square of the residuals, derived from 
equation (Al) and regarded as a function of the 
coefficients, is 


The components of the vectors are the n data points read 
from figure 2. We will assume that the x,- data is exact 
and that each 7 ,- point has an associated error A^,-. The 
norm of the y vector is defined by 

t,y] (A6) 


6 (Ci ,C3 — 



cik- \xi 




(A3) 


The condition that the sum of residuals be a minimum is 
expressed as 


~=0; 1=1 K (A4) 

Equation (A4) leads to a set of N linear equations 
which, with some matrix manipulations, can be solved 


The norm gives the length of a vector in the 
7i-dimensional space in which it is imbedded, and it is an 
invariant with respect to coordinate transformations. 

However, because of the errors Ay/ there is an 
uncertainty in the vector y 


N^ + ANy = 


LO/ + Ay, 




'/2 


=/(Ayi • • -Ay^) 


(A7) 
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In the last expression we concentrate on the norm as a 
function / of the errors Ay,-. Expanding / in a 
multivariable Taylor series in Ay,- and ignoring terms 
higher than first order result in 

+ AA^j, =/o + ( ^ )^ A;;, + . . . + )^Ay„ 

^y={-^)'Lyi^yi 

' i=i 

This can be written as a scalar product of two vectors 



where 



(A8) 


(A9) 


The relationship in equation (A8) is pictured 
geometrically in figure 16. 

The problem is to resolve the vector y into its 
components along a coordinate system which spans the 
/7-dimensional space. That is, we wish to find the 
subspace in which to expand y. We assume that all n 
vectors are not necessary. Furthermore, any component 
of y having a value less than the error in the norm of y 
will be considered ignorable. We proceed by defining 



il 

s 

• . x„)^. 

fli =Cl 

II 


02 =C3 

X3=ix\ . 

■ ■ 

fl3 =C5 

etc. 


etc. 


(AlO) 


where T denotes the transpose. Then, equation (A5) is 
rewritten as 


y= 

i 


(All) 


The set of vectors A",- span the entire space. 
Unfortunately, they are not orthogonal. Basically, this is 
the reason they are not convenient in finding the subspace 
of y. 

An orthogonal set of coordinates can be constructed 
using the Gram -Schmidt orthogonalization method. Let 
( IF] , 1 ^ 2 , . . .) denote the orthogonal set given by 


fVi=Xi 






W^=Xi 




X2*Wx X^»W2 
IFfIFi W2»W2 




W2 


etc. 


(A12) 


The data vector y can now be expanded in the orthogonal 
set as 


y=^bifVi (A13) 

t 

in which the expansion coefficients are given by 


bj=yWj 


(A14) 


and the right side is the scalar product of y and Wj. 
Finally, we impose the condition that 


Figure 16. - Geometrical reiationship between data ^ aaj 

vector 7 and error vector e. ^ 


(A15) 
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(A12), the highest order polynomial with which to apply 
the least-squares fit of equation (A4). 


and we reject all b’s which are less than ANy. 

Using this procedure we can define the subspace in 
which y lies. This also defines, because of equation 
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Appendix B 

Spectral Distribution of the Voltage Transfer Function and Normalization of the 

Spectral Power Density Function 


We assume that a normalized voltage transfer function 
(see eq. (Al)) has been found by least-squares fitting to a 
truncated polynomial series. It is also assumed that the 
truncation of the series is consistent with the inherent 
error of the data. 

Taking the Fourier Transform of equation (Al) results 


The integral of the spectral density over all frequencies 
can be carried out symbolically and the convolution 
integrals occurring in equation (B4) can be reduced to 
single integrations. Let 


Ey^r Y{f)df 

J — OD 


Y(f ) = ti ^2k - ‘] 


(Bl) Then 


where 


J — 00 


denoting 


Ey=t,C2k-S df 

*=1 


For notational clarity we change the variable of 
integration to 5i ; then the multiple integrations become 


X{f) = F\x(t)] (B3) 

We note that the Fourier transform of a product in the 
t-domain is a convolution integral of the transform in the 
/-domain. Hence, 

F{A2(r)] = r ;sr(/--5, )JT(5 i 1 ^X{f)®X{j) 

J — OD 


J — OD 

= J" ...r;<r(5i-52)^(^2-53)- 

J — OD J -00 


'XX(S2k-\)dS\ . . .dS2k- 1 


Fl^(/)]= X(f-Si)X(s^ -S 2 )X(s 2 )dsidS 2 


= X(f)®X(/)®X(f) 


-^(/)®3 


With this notation equation (Bl), the voltage spectral 
density, becomes 


The variables of integration are now changed using the 
following transformation: 


Sj =Mi +U2+ . . ■ -\-U2k-\ 


S2 =Ui +U2+ . . . +U2k-2 


S2k-l =“1 




I I 1 1 


I I HIM nil II 


The Jacobian of this transformation is unity; therefore, 
cfel . . . dsik - 1 = . . . dU2k - 1 (B9) 


This form assumes N Gaussian signals having 
amplitude a„ and centered about /„ . Each signal has a 
standard deviation of 1/V2X. By letting the index n range 
from -N X.O N (exclusive of zero) and allowing the 
convention 


and equation (B7) becomes 

J — OB 

= r ...rx{u2k-i)xiu2k-i)... 

J — oe J ^ Co 


=a„ 


(B14) 


f—n fnj 

equation (B13) can be rewritten more compactly as 


X(f) = 


£ One 

n=-V 




(B15) 


xX{ui)du2k-\. ■ -dux 



X(M2k-i)du2k-\ 


The prime on the summation will be taken to mean that 
rt =0 is excluded from the sum. The signal x(t) in the time 
domain corresponding to equation (B15) is found by 
taking the inverse Fourier transform F~^ of X(f) and is 


! 0 e 

XiU2k-2)dU2k~2- • • 

— OB 

{ Ob 

X{ui)dux 

OB 


-(r. ‘ 


(BIO) 


The last step, leading to equaton (BIO), results from 
the recognition that all the variables of integration M] . . . 
U 2 k-i are dummy variables and can be replaced by /. 
Equation (B6) can now be written as 


x(0=F-HX(f)] 

= r X(J)efi^^^df 

J —oo 

= 2( J) £ a„ cos (B16) 


Hence, the spectral distribution of equation (B15) results 
from applying the input signal of equation (B16) which is 
the application of N carriers for an effective duration, 
namely. 



(B17) 


(BID 

*=i 

with 


X(f)df (B12) 

J —00 

In order to proceed further with the spectral analysis, it 
is necessary to assume a particular form for X(f). We will 
assume that the input signal (double-sided, i.e., both 
positive and negative frequencies) is 


Here the effective duration is arbitrarily defined by the 
e-folding distance of the Gaussian time function. In fact, 
the duration is infinite. 

The response of the nonlinear system in the frequency 
domain is given by equation (B4) with equation (B15) as 
the input. In order to determine the general form of the 
2K-\ convolution integrals which are needed we will 
proceed with the calculation of the third-order 
convolution. Thus 

= Tt2 = -N n-^=-N 


X{f)=t<^n[ 

n = l *- 




(B13) 



XG(s„52,53)gj52rfS3 


(B18) 
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where 


G(Wi,W2,W3,W4,Ws) 


GiSi ,S2,S3) = 1 - 52 + (.S2 ~ S3 

+ (^3-//,)^ (B19) 

Using the following transformation: 

51 = «! +y^j ~^fn2 'I 

52 = U2+fni+fm ^ (B20) 

53=«3+/„3 J 

for which 
ds2dss =du2du3 
results in 


= + + + + (B25) 

and the seventh order becomes 

+ + + + (B26) 

It is clear how the higher orders can be transformed. The 
integrals in equation (B18), the third-order convolution, 
become 

r r ^-xaw3) ^,3 

J — 00 J — oo 


G(M],M2.«3) = («1 -«2>^ +(M2-M3)^ +W3 (B21) 


which can be written 


2 2 

G(«i,M2.W3) = 2(«3- 2“2) +2(“^~3"0 


Transforming variables again 


Wi =Ui 

2 

W2=U2~^Ux ^ 

1 

W’3=“3-2"2 


(B23) 


for which 



dw^ 


J — 00 

. '/i ^ 'A .«<» ,2 

=(s) (A) “'•) 

<kf (if 

= ('!) J_ (B27) 

VX/V3 


dw2dw3 =du2duj 


yields 


3 1 

G(wi,W2,W3) = 2w|-i- 2W^+ 3>vJ (B24) 

This procedure can be generalized for the higher order 
convolutions. The fifth order, for instance, is 
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The integrals for the fifth order become 

r r r r 

J — 00 J— oe J— 00 J— 00 

J — 00 J — 00 


/ 1 y/2 

^ 2 y/2 1 


( 4 N 


^3xi 

Uxi ' 

Isx/ 


= l(n%-(^/5).v? (B28) 

Vj vx/ 

And the multiple integrals from the seventh order become 


= t'- fV-~, 

V5 ni = -^ n^=-N 


xe 


j[/-U,+- +x,j)f 


(B32) 


t'--- L V--.7 

V7 ni^-N ny=-N 


xe 


-^[/-On,+ - +/.7)f 


(B33) 


And in general the /sTth convolution is 


Xe 


-/-(/■«,+ ■•■ +//ij] 


(B34) 



Finally we list the convolutions of various order needed 
for the evaluation of equation (B4). In these expressions 
we have resubstituted variables according to 
W] =Mi =/-/„! +/rt 2 + • • • +fnfc where K is the Ath 
convolution. Thus 


The evaluation of these multiple sums is the subject of 
appendix C. As a final note we observe that the 
integration over all frequencies will yield the same results 
as in equation (BIO). 

The final topic of this appendix is the normalization of 
the spectral power density function. By definition the 
square of equation (B4) is the spectral power density; that 
is. 


fC K 

Atj = 1 * 2=1 


XV)- 


«1= -iV 


xm^^-j^l t' t' t'w 

V3 ^ /f| = — nj = — A/ h2=—N 


"3 


xe 


2 f 


X{f)®(?-ki - 1)A^(/)®(2*2 - 1) (B35) 

(B30) 

Referring to equation (Bll) note that at saturation Ex 
and Ey are both unity because the voltage transfer curve 
(fig. 2) has been normalized. Therefore, 

f;c2*:-i=l (B36) 

t=i 

(B31) The power gain curve of figure 1 must be derivable 
from the spectral power density by integrating over all 
frequencies; thus. 
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And, 


K K 

r n{f)df= £ E C2k,-ic2k,-i 

*1 = 1 *2=1 

r (B37) 


Because the voltage transfer curve was derived by taking 
the square root of the power curve, equation (B37) must 
be expressible as a perfect square having the following 
form: 


K K 

r Y^{f)df= 2 i; C2k,-ic2k,-is 

*1 = 1 k2=\ 

K 

E 

*1=1 


-2{A:,+Ar2-1) 


*1 = 1 *2 = 1 

= ( E ^ 2 *- lS 2 *:-iy 


,2/t2-l 


(B43) 


At saturation d = 1 and, therefore, the right side is equal 
to one. 

The convolutions in equation (B37) can be expressed in 
terms of the convolution of saturation by means of a 
scale factor S; that is, 

^(/)®K=SKX,at(/)®K (B39) 


This is most easily seen by recalling the origination of 
these terms 


At saturation, S= 1 and equation (B43) because of 
equation (B36) becomes 

(“ Y^{f)df=\ (B44) 

Furthermore, the integral oi Y^{f) between two finite 
frequencies /a and ft, is 

f Y^{f)df= f t ^2*1-1 

'‘fa *1 = 1 *2=1 

XC2it2-iS2^*'^*^"‘^7’*i.*2 (B45) 


^(/)®K = F[xfi(t)] = S^/^ATfat (0] (B40) 


With these scale factors equation (B37) can be written as 

r Y^Wf=t E 

*1 = 1 *2= * 

xf (B41) 

J — oo 


In order that equation (B41) be reducible to the perfect 
square form of equaton (B38), it is clear that each term of 
the expression for Y^(f)hG normalized by the integral of 
the term at saturation. Hence, 


where 


^*l,*2 “ 


r ^ Q[2 Ati - 1) Q[2A:2— 1)^^ 

I sat sat 



J — OD 


(B46) 


Equation (B46) expresses the interesting result that the 
entire nonlinear problem is reduced to an algebraic 
relationship in terms of its characteristics at the 
saturation point. This allows a great simplification in 
computer coding and a savings in computer time. 


J^(/)= E E <^2k,-lC2k3-I 

* 1=1 * 2=1 

X [aT(/)®(2^. - 1) AT(/)®(2<^2 - 1)/ 

J (/)®(2*> - - *>^/] (B42) 
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Appendix C 


Evaluation of the Multiple Sums Occurring in the Output 


Spectral Density Function 


The first step in the evaluation of the multiple sums 
occurring in equations (B31) to (B34) is to simplify the 
notation to 




fl ~ =*=(//! I 


(C6) 


II 

® 

L' 

E' £'^«1+/«2+/«3) 

/i2== -N ny= --/V 

(Cl) 

Defining 


X(f)®^ = 

^ # 

ni = -A' 

... E' ■■■ +/.s) 

/7j= -A/ 

(C2) 

Wi) = F(/i)-H/X-/,)^ 

H(f2)=m)+n-f2U 

(C7) 

X(f)®^ = 

^ I 

L’ 

ni = ~N 

... E' HTni+ ... +/n/t) 

n^=-N 

(C3) 

Equation (C5) can be simplified to 


Note that equation (Cl) crm be written as sums over 
positive integers alone. Thus, the sum over nj becomes 

_ N N N 

XV)®=. E E E 

rtj = l ^2=1 «3=l 

~^fn2 “^.^3) +^n2 ~fn^ 

(C8) 


_ N , N N 

n^=— N ri2=-N ^3 = ! 

~^fn2 '^fn2 


(C4) 


Resolving the sums over n \ and ri 2 in the same way yields 


nj=l /i2=l n3 = l 

~^fn2 ~fn\ ~^fn2 '^fn^ 

<t 1 ~fn 2 ~fn j ~fri 2 '^fn^ 

■h f^n\ ■f/n2 ~fn^ ~fn\ ~fn^ 

+ 1 1 ~fri2 ~fn-^ ~fn j ~fri2 ~/n 3 )] (^ 5 ) 

Because n\, ri 2 , are dummy indices whose only 
purpose is to provide distinctness, they can be 
interchanged at will. Hence, equation (C5) contains only 
two combinations of frequencies: 


The coefficient 3 in the last term arises from counting the 
number of ways one negative sign can be distributed over 
the three frequencies given that permutations of the two 
positive signs are not distinct; that is. 



(C9) 


Actually the coefficient of the first arises in the same way 



(CIO) 


That is, it is the number of ways to distribute zero 
negative signs over the three frequencies divided by the 
indistinct number of permutations of the positive sign. 

In the case of the fifth-order sums, the same procedure 
can be applied and gives the following result: 
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(C16) 


It is clear from these two cases how the higher order sums 
can be rearranged. 

The purpose of this exercise is to reduce the number of 
calculations that must be made. The total number of 
frequency computations in equations (Cl) and (C2) is 
(2A/)^ and (2N)^. Of these, most are redundant. The 
rearrangements given in equations (C8) and (Cll) reduce 
these to 2N^ and 3N^ a factor of 1/4 and 3/32, 
respectively. For the seventh-order sums, the reduction 
factor is 1/32; and for the ninth order, only 5/128 of the 
total calculations are required. However, the remaining 
calculations are still many and further reductions sure 
needed. These can be made through the narrow band 
approximation which reduces the calculations from 
(2N)^ to NJ<. 


fn\~^ fn2 ^n2 ^ 

Obviously, the combination of integers {k„^ -i- k „2 - k„^) is 
also an integer. The same is true for orders higher than 
third. 

If the structure of input frequencies is complicated 
and/or the order of intermodulation products is high, it 
is convenient to use a computer to calculate the 
frequencies in equation (Cl 4). A computer program has 
been written in which the various orders are computed 
recursively. The program coding will not be given here, 
but its computational procedures will be outlined. 

First, a second-order set of integers is computed; 
namely. 


Let 

/=/o + " (C12) 

where 

-W<p:<]V (C13) 

where 2 IF is the bandwidth. In the narrow band 
approximation the computations are limited to the band 
of frequencies centered on /o . Therefore, the first term in 
equation (C8) can be ignored because its argument is 
centered on 3 /q. Likewise, the first and second terms of 
the fifth-order convolution in equation (Cll) can be 
ignored because they are centered on 5/o and 3/o- For 
each of the convolutions the terms that are retained have 
the following arguments: 

Third order: +f „2 


-ic 


ni =1, 


«2 = 1. 


N 


,N 


(C17) 


Associated with these numbers is an amplitude factor 
given by 



*rt2 


(C18) 


The range of i in both of these equations is clearly 
between 1 and N 2 ■ 

The third-order set of integers and the associated 
amplitude factor can now be computed 


Fifth order: /„, +f „2 +fn-i -fn^ ~fni 


Seventh order: /„, +/„2 +X13 -f„^ -/„g (C14) 

etc. 


K^^)=k„,+Kp-^ 


ni =1, . . . , N 

/=1 


(C19) 


In order to compute the resultant frequencies given in 
equation (C14), we let 

A=/o+V (C15) 


If the calculations in equation (Cl 9) were carried out as 
indicated, then there would be integers and amplitude 
factors. However, not all of the integers are unique. 
For each redundancy the amplitude factors are added and 
become a weighting function. Hence, 


where w is an appropriate increment and A:„. is an integer. 
Such a conversion of frequencies to integers can always 
be accomplished with arbitrary accuracy. The third-order 
frequencies (from eq. (Cl 4)) can then be written 


aS' = E Ai^y 

m‘ 


(C20) 


21 



The range of m is from 1 to M^) in which M3) is to (B34), (C12), (CM), and (C15) and considering only 
considerably less than M . the positive frequencies result in the coherent terms being 

The fifth-order terms are now computed from 


k^^)=kP+kP ] /=1 M 

I y = l, ...,M3) 

m ^ 


where 

m = l, , M3) < <M 


n= 1 

(C2i) 






-(X/3X»'-^„>v)2 


Again, in the computations of equation (C21) the 
accounting for redundancy is indicated by the sum over 
w and results in M^) values. 

The seventh-order terms are 


=/irP +kP 


< = 1, . 

. . ,M 


► 

7=1. ■ 

. . , A(3) 

m ^ 




where 





(C22) 



a5e-(X/5)(v-Ar„H-)2 




(\/7K^-k„w)2 





(C24) 

(C25) 

(C26) 

(C27) 

(C28) 


m = l, . . . , M'^) < <A^ 


And, finally, the ninth -order terms are 


And the general term where K is an integer is 


aS'‘E'>FW’ ' 

m J 


/=l, . . . ,M 

y=l, . . . ,M’) 


where 

M9)<<n9 


(C23) 


X^(J)®{2K+l)= 

\/2AT1 

X ^ O p/f + i)g - H» - k„wm2K+ 1) ^^^29) 


The program is limited to the ninth order but it can be 
easily extended to higher orders. 

The purpose of these calculations is to compute the 
amount of signal distortion that is introduced by the 
nonlinearity. But not all of the nonlinear terms in 
equation (B4) can be considered as noise. The higher 
terms do contribute to a coherent amplification of the 
input signal. If we assume that the information will be 
carried on the amplitude factors, then crosstalk tmd other 
distortions can be attributed to terms involving cross 
products of amplitudes. Therefore, using equations (B30) 


Finally, the complete expressions for the convolutions 
including the coherent terms are as follows: 


^ 1 


Xif) 


®5 _ 


1,2 /v(J) 


VJVX 


m= 1 


(C31) 
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-f- 

V7^X 


3 /v(7) 


) 


(7)« 


-X(;-A<^>h.)V7 


2 


4 ^/») 


(D 2^ 


(9)f 


-X((--/i^’'w)V9 


And the general K\h term is 

(C32) 


'IW+I 



(C33) 


M^+i) 

X ^ y^(2K+\)^-Hv-f^*»w)Y/(2K+\) 
m - 1 


(C34) 
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